DUDLEY .  LiE.sARY 

NAVAL  POSTGRADUATE  SCHOOL 
MONTERcY,  CALIF.   93940 


An  Algebraic  Propeller  Model 
Using  Eulerian  Turbomachine  Theory 

by 

Dennis  Ray  Milken 

B.S.E.E. ,  Purdue  University 
<  1974) 

SUBMITTED  TO  THE  DEPARTMEI^  OF  OCEAN 
ENGINEERING  IN  PARTIAL  FULFILLMENT  OF  THE 

DEGREES  OF 

OCEAN  ENGINEER 

AND 

MASTER  OF  SCIENCE  IN 
MECHANICAL  ENGINEERING 

at  the 

MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 

MAY  1982 

©   Dennis  Ray  Nil  ken   1982 

The  author  grants  to  M.I.T.  permission  to  reproduce 
and  to  distribute  copies  o-f  this  thesis  document  in 
whole  or  in  part. 

/ 


DUDLEV  ,  :     .  .  -  '.Y 

NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY,  CALIF.   93940 

An  Algebraic  Propeller  Model 

Using  Eulerian  Turbc»nachine  Theory 

N  by 

Dennis  Ray  Milken 

Submitted  to  the  Department  o-f  Ocean  Engineering  on  May  7,  1982 
in  partial  -fulfillment  o-f  the  requirements  -for  the  degrees  o-f 
Ocean  Engineer  and  Master  o-f  Science  in  Mechanical  Engine- 
ering. 

ABSTRACT 

Utilizing  -four -quadrant  experimental  data  determined  by  Minio- 
vich  converted  to  the  modified  advance  and  angl e-of-advance 
coefficient  form,  complete  Karman-Knapp  Circle  Diagrams  are 
plotted  for  a  number  of  3  bl aded  propellers.  Reducing  the 
family  of  curves  to  a  rational  basis,  algebraic  equations  are 
derived  to  model  the  performance  of  a  three  bladed  propeller 
through  all  four  quadrants  using  Eulerian  turbomachine  theory. 
The  result  is  a  set  of  equations  uihich  can  be  used  to  generate 
complete  propeller  thrust  and  torque  characteristics  given  only 
incomplete  performance  data  suitable  for  use  in  steady  state 
design  or  dynamic  simulation. 
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CHAPTER  1      INTRODUCTION 

1 . 1  Primary  motivation  for  this  kind  of  modelino  method 

For  a  naval  architect  or  anyone  interested  in  the  modeling  of 
propellers  as  a  part  of  the  design  procedure  leading  to  an 
integrated  ship  design,  the  problem  long  persisted  where  one 
is  required  to  estimate  the  performance  of  a  propeller  for 
ship  modeling  purposes  before  the  design  has  progressed  to  the 
point  where  one  is  reasonably  confident  as  to  what  the  final 
desired  propeller  should  look  like.  Even  more  frustrating  is 
the  problem  of  having  somehow  determined  the  geometry  of  the 
propeller  but  not  having  sufficient  experimental  data  to  allow 
a  comprehensive  computer  simulation. 

Even  where  it  has  existed,  considerable  effort  has  been  expen- 
ded to  convert  the  data  into  useable  form  for  ccwnputer  simu- 
lations since  the  practice  of  tabulating  the  results  as  func- 
tions of  an  advance  coefficient  has  difficulty  in  the  vicinity 
of  zero  ship  or  shaft  speed. 

The  problems  associated  with  modeling  a  machine  such  as  a 
propeller  are  formidable  since  a  complete  model  must  make 
allowances  for  complicated  fluid  motions  which  are  only  just 
beginning  to  be  understood.  In  particular,  it  is  still  quite 
difficult  to  pass  completely  from  propeller  geometry  to  the 
detailed  lift  and  drag  forces  required  for  specific  propeller 
characteristics. 

However,  there  now  exists  a  considerable  amount  of  documen- 
tation  for  modeling  machines  y^ery   much  like  a  propeller.   For 


example,  industrial  axial  turbines  and  pumps  have  had  the 
attention  of  modelers  and  those  who  perform  simulations  -for 
decades.  Taken  together,  propellers,  turbines  and  pumps  are  in 
fact,  members  of  a  general  class  of  machines  known  as  turbo- 
machines.  Models  for  many  of  the  large  industrial  machines 
have  been  quite  successful  in  predicting  machine  performance 
characteristics  and  are  widely  used  in  control  system  synth- 
esis and  design.  A  general  model  for  propellers  which  could 
incorporate  that  data  base  would  have  a  firmer  basis  in  exper- 
ience and  first  principles. 

1.2  Some  current  propeller  model i no  methods 

At  present  there  are  two  common  methods  used  by  designers  to 
obtain  propeller  data  for  computer  simulations.  The  first 
involves  experimental  determination  of  the  characteristics  of 
a  physical  model  of  the  full  size  propeller  followed  by  con- 
version of  the  data  into  the  familar  thrust  and  torque  coeffi- 
cient form.  Such  coefficients  are  then  used  in  a  table  look-up 
scheme  within  the  simulation  program  to  obtain  the  necessary 
thrust  and  torque  values.  There  are  advantages  to  this  kind  of 
method.  First,  the  determination  of  the  propeller  characteris- 
tics through  experimentation  (scaling  considerations  aside) 
ensure  that  the  data  used  in  the  simulation  will  be  as  accu- 
rate as  possible.  Secondly,  the  use  of  a  model  of  the  intended 
propeller  acts  as  a  test  of  the  geometry  and  reveals  at  an 
early  date  any  serious  conceptual  problems. 
The  disadvantages  of  such  a  method  are  fairly  obvious.  To  test 


a  model  o-F  a  propeller  at  an  early  stage  in  a  design,  you  must 
have  a  well  defined  idea  as  to  the  range  of  operation  of  the 
propeller  and  its  desired  geometry.  To  allow  for  a  comprehen- 
sive simulation,  you  must  test  and  record  the  model  data  in 
all  four  quadrants  of  operation.  Since  cost  is  always  a  signi- 
ficant consideration,  the  machining  of  a  propeller  model  and 
the  set-up  and  operational  costs  of  the  tests  represent  a 
sizeable  investment  in  the  design.  As  alluded  to  earlier,  raw 
data  is  not  suitable  for  computer  simulations  due  to  the 
problem  of  "singularities'  of  the  advance  coefficient  at  zero 
ship  speed  and  zero  shaft  rpm.  The  experimentally  determined 
data  must  be  converted  to  a  form  that  avoids  such  difficul- 
ties. 

The  second  method  of  modeling  propeller  characteristics  invol- 
ves complex  and  relatively  recent  work  in  propeller  hydrody- 
namics. The  computer  models  of  Professors  Kerwin  [13  and 
Cummings  C23  are  quite  thorough  and  provide  some  revealing 
results.  As  a  "complete"  model,  these  computer  programs  are 
undoubtedly  at  the  vangard  in  propeller  theory.  With  a  firm 
basis  in  hydrodynamic  theory,  the  continual  evolution  of  these 
programs  will  ultimately  achieve  very  accurate  results.  Unfor- 
tunately, these  programs  are  imnensely  long  and  require  consi- 
derable amounts  of  CPU  time.  Because  of  the  attendant  expense, 
they  are  not  yet  practical  for  incorporation  directly  into  a 
ship  simulation  program. 
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1 .3  The  Objective  o-f  the  Modeling  Equations 

The  objective  o-f  this  thesis  will  be  to  derive  a  set  o-f  rela- 
tively simple  algebraic  equations  that  will  serve  as  a  tool 
for  naval  architects  in  the  ship  design  process.  By  being  able 
to  predict  accurately  propeller  per-formance  in  all  -four  quad- 
rants o-f  operation  based  on  Knowledge  obtained  from  limited 
experiments  and  in  a  form  suitable  for  computer  simulation, 
the  naval  architect  can  save  considerable  time  and  expense  in 
the  simulation  portion  of  the  design.  Experimental  determina- 
tion of  the  points  of  minimum  dissipation  would  be  sufficient 
to  allow  for  a  comprehensive  simulation  for  forward  or  crash- 
back  maneuvers.  Such  a  set  of  equations  may  not  eliminate  the 
need  for  actual  experiment  in  the  final  stages  of  the  design, 
however  they  would  greatly  simplify  the  initial  requirements. 
Additionally,  where  experimental  data  exists  for  particular 
propellers  in  only  the  first  quadrant,  the  equations  could  be 
used  to  reconstruct  the  remaining  three  quadrants  with  good 
accuracy. 


CHAPTER  2      DERIVATION  o-f  the  ALGEBRAIC  MODEL 

2.  1  Fundamentals  of  Eulerian  turbcwnachines 

The  practice  o-f  presenting  open  water  propeller  per-formance 
data  as  nondintensional  thrust  <Kt)  and  torque  (Kq)  coe-f-fi- 
cients  versus  an  absolute  advance  coefficient  (lambda)  or 
simply  advance  coefficient  <J)  C3],C4]  incorporates  assump- 
tions which  inherently  restrict  the  machines  of  interest  to 
Eulerian  turbcxnachines.  That  is,  by  considering  only  high 
Reynolds  numbers  the  performance  of  the  machine  is  governed  by 
inertia  and  pressure  forces.  Dimensional  analysis  on  the  ex- 
pression : 

(Thrust)     F=f<'v»p,D,N,p) 
yiel ds 

Kt=  F/PD4N2=f  (^v'p/ND)  =f<J)  <2.1) 

Similarly  for  torque, 

Kq=  T/pD5N2=g<v^p/ND)  =g<J)  <2.2) 

where  Vp  is  the  speed  of  advance  of  the  propeller  through  the 
water  and  D  any  characteristic  propeller  dimension  (usually 
diameter).  This  single  parameter  representation  of  propeller 
performance  yields  the  powerful  result  that  the  performance 
data  of  geometrically  similar  propellers  will,  when  plotted 
together,  generate  to  good  approximation  a  single  curve. 
The  same  assumptions  would  ,  from  a  fluid  mechanics  stand- 
point, give  Euler's  equation  of  inviscid  motion.  That  is, 

pv3eGRAD(v)  =  -6RAD(P)  (2.3) 
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Or,  in  a  more  general  -form  to  incorporate  unsteady  -flows, 

PD<v)/Dt    =    -6RAD<P>  <2.4) 

where  D/Dt  is  the  substantive  derivative. 

In  conjunction  with  equation  (2.4),  we  express  mass  conserva- 
tion by: 

D<P)/Dt  +  PGRADJ^v  =0  (2.5) 

where  P  is  independent  o-f  position  and  v  is  the  vector  -fluid 
velocity.  Re-ference  C5]  chapter  3  and  Re-ference  161  are  use-ful 
for  a  more  detailed  analysis. 

A  mach  ine  satis-fying  equat  i  ons  <  2.  4)  and  (2.5)  will,  as  ver  i- 
-fied  experimentally,  demonstrate  what  is  re-ferred  to  as 
■Eulerian  similitude".  That  is,  geometrically  similar  machines 
will  produce  similar  per-formance  characteristics  when  present- 
ed in  a  properly  chosen  nondimensi  onal  -form.  Casanady  C73 
re-fers  to  this  as  a  'per-formance  law'. 

Based  on  the  discussion  above,  the  determination  o-f  the  per- 
■formance  characteristics  o-f  a  large  turbcwnachine  is  usually 
done  with  small  scale  models.  The  -full  scale  results  are 
derived  -from  the  model  results  by  recognizing  that  at  corres- 
ponding points  on  a  characteristic  curve,  all  nondimensi onal 
parameter  ccmbinations  are  also  equal. 
However,   there  are  problems  with  scale  e-ffects  not  evident  in 
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such  a  single  parameter  representation.  In  addition  to  geo- 
metric dissimilarities,  such  as  roughness  and  clearances,  hy- 
drodynamic  differences  are  also  present.  These  include  Rey- 
nold's and  Froude  number  variations  as  well  as  cavitation  and 
compressibility  number  differences.  Attempting  to  incorporate 
some  or  all  of  these  variations  into  a  performance  law  would 
require  that  the  law  be  multidimensional.  As  an  example,  one 
could  attempt  to  formulate  a  nondimensional  representation  for 
thrust  and  torque  characteristics  by: 

Kt=f<J,Re)      Kq=f<J,Re)  . 

Here,  the  equivalence  of  the  advance  coefficient  in  conjun- 
ction with  the  Reynolds's  number  would  lead  to  similitude. 
Similarly,  one  could  also  form  yet  another  expression  for 
thrust    and    torque 

Kt=f<J,P)  Kq=f<J,P> 

where  P  represents  a  nondimensional  ratio  for  propeller  pitch. 
Such  a  formulation  could  be  applicable  to  Kaplan  turbines  or 
controllable  pitch  propellers  where  the  blade  pitch  is  varied 
mechanical  1 y . 

Multi-dimensional  performance  laws  unfortunately  lead  to  fami- 
lies of  curves,  a  problem  which  single  parameter  laws  avoid. 
As  an  example,  Bovet  C8]  presents  efficiency  versus  output  at 
a  constant  head  for  a  Kaplan  turbine  for  various  blade  set- 
tings. An  analytical  treatment  for  a  multi-dimensional  per- 
formance law  in  terms  of  generating  numerical  data  would 
entail    a  matrix    of    possible   states  with    dimension    equal     to    the 

12 


number  o-f  nondimensional  characteristic  parameters.  Obviously, 
dimensional  analysis  techniques  used  to  generate  characteris- 
tic nondimensional  performance  parameters  are  limited  in  their 
applicability  in  performance  simulation.  To  allow  incorpora- 
tion o-f  fluid  or  geometrical  affects  into  the  analysis  Mould 
require  another  approach.  That  approach  should  produce  a  for- 
mulation which  retains  the  advantages  of  similitude  yet  allows 
for  some  fluid  or  geometrical  effects. 

2.2  Ideal,  canonical  and  perfect  ship  propeller  equations 

Figure  2.1  shows  a  propeller  blade  cross  section  at  some 
arbitrary  radial  distance  from  the  propeller  hub  under  the 
condition  of  zero  nominal  speed  of  advance.  Figure  2.1a  in- 
dicates the  tangential  velocity  of  the  blade  element  (Vr)  and 
fluid  velocity  (Vf ) .  Figure  2.1b  shows  the  lift  (L)  and  drag 
<D)  forces  with  respect  to  an  effective  angle  of  attack  a.  The 
magnitude  of  the  effective  angle  of  attack  is  a  function  of 
the  circulation  which  causes  an  induced  velocity  [33. 
If  we  define  a  torque  exerted  by  the  element  as  Fr  then, 

Fr=  LSin<oc)  +  DCos<a)  <2.6) 

Further,  defining  Ff  as  a  force  perpendicular  to  Fr  in  the 
direction  of  fluid  motion  (fluid  pressure  term), 

Ff=  LCosCo)  -  DSinCo)  (2.7) 
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<a) 


Fi  Qure   2 . 1 
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Since  L  and  D  are  expressed  as: 

L=   Cl<i4pAvXv)     ;    D=   Cd<J4pAvXv)  , 

<2.6)    and    <2.7)    can    be   rewritten: 

Fr=   A^pAvCClX^f    +    CdJfVr)  (2.8) 

F-f=   i4PAv(CdXVr    -   ClXV-f)  <2.9) 

In    per    unit    (see  Appendix   A),    (2.8)    and    (2.9)    become: 

m=  v(ClXq   +    CdJfn)  (2.10) 

p=  v(ClXn    -    CdJCq)  (2.11) 

Or  using  the  canonical  notation, 

m=   RgJ^q   +    RdSCn  (2.  12) 

p=   Rg3«n    -   RdJfq  (2.  13) 

where 

Rg=rq3ev  rq=rq(P)  v=SQRT(nXn  +  q3Cq) 

Rd=rdJev  rd=rdO)  ^=ARCTAN(q/n) 

The  correspondence  with  the  bond  graph  representation  shown  in 
Appendix  B  is  now  readily  apparent. 

The  mechanical  power  can  be  expressed  as  Fr^Kfr   while  the  -fluid 
power  can  be  expressed  as  F-fX'v'-f.   Thus,  dissipation  losses  can 
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be   -Found    to   be 


Pd=   Fr3?Vr    -   FfXV-f  (2.14) 


In   per    unit,    (2.14)    becomes 


Pd=  m^n    -   p3Sq  <2.  15) 

Substituting  (2.12)  and  <2.13)  yields 

Pd=  RdMvS^v  <2.  16) 

Hence,  as  would  be  expected,  an  ideal  ship  propeller  would 
have  no  power  lost  in  dissipation  i-f  it  could  generate  lift 
without  drag.  The  equations  reduce  to  the  rather  trivial 
expressions 

m  =  Rg3?q  <2.  17a) 

p  =  RgXn  <2. 17b) 

with  Rg  representing  an  energy  trans-formation  modulus.  Rd  can 
be  considered  as  a  general  loss  term  that  includes  a  combin- 
ation o-f  a  number  o-f  loss  mechanisms  such  as  -form  drag  and 
separation  losses.  Writing  Rd  in  a  more  general  -form  typical 
o-f  such  losses  such  that  at  the  point  o-f  minimum  dissipation 
loss,  Rd  is  identically  zero  [93,  we  can  de-fine  a  PERFECT 
propel  1 er . 
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Rd  =  W<P)3f<<n-q)5«<n-q)/v)  <2.18) 

Using  the  Kutta  condition  to  -formulate  a  general  expression 
-for  the  modulus  Rg  would  suggest 

Rg  =  rg  -»^  XO)  3«(  <nXn) -<qXq)  )/v  <2.19) 

di-f-fering  -from  the  -form  proposed  by  Paynter  by  the  term  rg 
which  could  represent  some  linear  -function  in  n.  Substituting 
<2.18)  and  (2.19)  into  (2.12>  and  (2.13)  and  recognizing  the 
arbitrariness  in  WO)  and  XO)  ,  we  can  arrive  at  the  -following 
expressions: 

m  =   rgSSq    *    Kl<P)5Cv<n-q)  <2.20) 

p  =   rgXn    -^    K20)  5^v<n-q)  (2.21) 

where  Kl  and  K2  now  incorporate  the  individual  character  o-f  a 
given  propeller.  The  equations  above  -for  a  PERFECT  propeller 
(zero  minimum  dissipation  loss)  can  be  au^nented  to  represent 
a  PRACTICAL  propeller  by  including  an  additional  term  to 
re-flect    the   actual    non-zero  minimum   dissipation    loss. 

m  =   rg^q   +    KKP)  3fv<n-q)    *    R13^n  (2.22) 

p  =   rgXn    -»•    K2 ( P)  3C v ( n -q)    -   R23^q  (2.23) 

Rl      and   R2  would   be   -functions   o-f    the   advance   coe-f -f icient    J   -for 
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a  given  propeller,  but  in  a  simpler  model  could  be  reduced  to 
constants.  As  will  be  seen  later,  the  values  o-f  Rl  and  R2  have 
a  direct  effect  on  the  shape  o-f  the  4-quadrant  curves  o-f  a 
given  machine  and  in  particular,  the  slope  o-f  the  zero  torque 
and  thrust  radial  s.  KIO)  and  K20)  -for  an  idealized  propeller 
would  be  identical.  The  degree  o-f  agreement  between  Ki  and  K2 
based  on  experimental  data  would  be  a  measure  o-f  the  validity 
o-f  the  model  -for  a  given  rg. 

2.3  Determination  o-f  the  algebraic  model  ina  equations 

Using  equations  <2.22)  and  <2,23)  as  a  model  -for  determining 
propeller  characteristics  requires  a  determination  o-f  the  K<P) 
term  as  well  as  Rl,  R2  and  rg  based  on  experimental  data.  The 
data  produced  by  Miniovich  CIO]  and  converted  to  a  modi-fied 
advance  coe-f-f icient  -form  Cll]  was  used  in  an  abbreviated 
methodology  in  order  to  -formulate  a  general  model  that  could 
be  expanded  upon  in  a  more  detailed  analysis.  See  also  Refer- 
ence C 12] .  Such  an  analysis  would  include  an  investigation  of 
those  factors  which  collectively  yield  the  individual  charac- 
ter of  a  given  machine  yet  have  defineable  specific  influen- 
ces, the  details  of  which  are  beyond  the  scope  of  this  thesis. 
In  Reference  Cll],  experimental  results  for  thrust  and  torque 
are  tabulated  for  eighteen  three  bladed  propellers  with  vary- 
ing expanded  area  ratios  and  pitch-to-diameter  ratios.  The 
propeller  designations  are  those  used  by  Baker  and  Patterson. 
Table  2.1  is  provided  as  a  quick  reference  for  EAR  and  P/D 
values  for  each  propeller. 
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Figures  B3  and  B4  in  Appendix  B  are  plots  o-f  rauj  data  for  two 
series  o-f  propellers.  To  be  useful,  these  curves  must  be 
rationalized  to  a  consistent  basis.  To  do  so,  the  points  o-f 
minimum  dissipation  Mere  -found  for  each  of  the  propellers  (See 
Appendix  A> .  Those  points  were  then  used  to  scale  figures  B3 
and  B4  such  that  the  point  of  minimum  dissipation  in  the 
region  of  normal  pump  operation  (quadrant  !>  was  fixed  at 
(l,i>  and  in  the  region  of  normal  turbine  operation  (quadrant 
III)  at  (-1,-1)  in  the  q  vs  n  plane.  The  mapping  of  the 
data  to  those  points  of  minimum  dissipation  is  done  over  the 
entire  q  vs  n  plane  recognizing  the  'mirror  image''  pattern  of 
the  curves.  That  is,  the  point  of  minimum  dissipation  in  the 
first  quadrant  for  positive  torque  and  thrust  is  reflected 
into  the  third  quadrant  on  the  negative  torque  and  thrust 
curves.  Hence  the  scaling  is  reflected  through  the  origin 
rather  than  simply  across  the  zero  torque  and  thrust  radial s. 
Figures  2.2  and  2.3  are  plots  of  the  thrust  and  torque  curves 
respectively  for  three  propellers  selected  for  their  varying 
expanded  area  and  pitch-to-diameter  ratios  and  are  representa- 
tive of  the  form  of  the  4-quadrant  curves  for  all  of  the 
propellers.  With  the  data  plotted  in  this  form,  the  similarity 
of  the  torque  and  thrust  curves  versus  the  advance  angle  is 
striking.  In  addition,  the  general  form  of  the  curves  is  quite 
consistent  in  the  regularity  of  the  specific  perturbations 
which  influence  all  of  the  propellers  in  the  same  manner. 
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Figure    2.2 
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Figure    2.3 
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Prop  No. 12  3           4  5  6 

EAR            0.5      0.5  0.5  0.5  0.5  0.5 

P/D            0.6      0.8  1.0  1.2  1.4  1.6 

Prop  No. 7           8  9  10  11  12 

EAR           0.8      0.8  0.8  0.8  0.8  0.8 

P/D           0.6      0.8  1.0  1.2  1.4  1.6 

Prop   No. 13          14  15  16  17  18 

EAR            1.1       1.1  1.1  1.1  1.1  1.1 

P/D           0.6      0.8  1.0  1.2  1.4  1.6 

Table  2. 1 


Be-fore  continuing,  I  must  point  out  that  there  is  an 
approximation  which  was  used  which  requires  an  explanation.  As 
Figure  B2  in  Appendix  8  demonstrates,  the  lines  o-f  zero 
thrust  (head)  in  the  -first  and  second  quadrants  are  obviously 
not  180  degrees  apart.  The  di-f-ference  between  the  angle 
measured  clockwise  -from  the  first  quadrant  zero  thrust  radial 
and  the  second  minus  180  degrees  I  call  the  skew  angle.  The 
skew  angle  -for  high  spec i -fie  speed  axial  -flow  pumps  is 
generally  quite  small.  See  re-ference  [133.  The  skew  angle 
computed  -for  the  propellers  used  in  this  analysis  were  all 
less  than  3  degrees  and  were  generally  less  than  1  degree.  The 
amount  o-f  skew  angle  is  a  -function  o-f  the  magnitude  o-f  losses 
o-f  the  machine.  Since  the  skew  angle  is  quite  small,  the 
curves  were  approximated  to  a  PERFECT  propeller  by  dividing 
the  skew  angle  between  the  two  radial  s  -found  -from  the   minimum 
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dissipation  calculations.  The  scaling  values  for  q  and  n  were 
then  recalculated  using  these  adjusted  angles.  As  a  result, 
the  -first  and  third  quadrant  asymptotes  were  180  degrees 
apart.  The  approximation  was  -felt  Justified  since  the  calcu- 
lated minimum  dissipation  angles  were  obtained  with  a  confi- 
dence of  plus  or  minus  2  degrees  due  to  the  required  extrapol- 
ations in  the  corresponding  portion  of  the  experimental  data. 
Such  an  approximation  also  simplifies  the  derivation  of  the 
modeling  equations  considerably.  Table  2.2  lists  the 
calculated  minimum  dissipation  angles  and  the  approximate 
angles  in  both  quadrants  for  six  of  the  eighteen  propellers 
along  with  the  calculated  and  approximated  magnitudes  of  the 
thrust  and  torque  radial s  at  those  angles.  The  approximate 
values  are  shown  in  parenthesis. 

Prop.  No.      Min.  Dissipation  Radial  Length 

angles  thrust        torque 

1  35.2  <   34.6)       11.52(8.97)   15.39(14.88) 
-146.0  (-145.4)       3.92(3.80)   25.83(22.44) 

2  39.5  (   38.5)       5.64(5.02)   13.54(11.32) 
-142.5  (-141.5)       4.98(4.59)   25.65(23.88) 

3  44.5  (   44.3)       4.95(4.82)   10.54(  9.96) 
-136.0  (-135.8)       5.02(4.90)   19.99(18.68) 

5  53.5  (   53.0)       5.06(4.78)    9.27(  8.85) 

-127.5  (-127.0)       6.88(6.31)   26.99(20.92) 

15  45.6  (   45.6)       6.99(6.90)   12.10(12.04) 

-134.5  (-134.5)       4.82(4.79)   20.36(19.61) 

18  58.0  (   57.8)       6.99(6.59)   11.96(11.47) 

-122.5  (-122.3)       7.48(7.07)   20.99(18.76) 

Table  2.2 
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Using  equations  <2.20)  and  (2.21)  for  the  PERFECT  propeller 
and  setting  rg  equal  to  n,  we  can  plot  the  K<P)  -functions 
required  to  -fit  the  experimental  data  by  setting  m  and  p  to  a 
constant  (1)  and  solving  -for  Kl  and  K2  at  each  data  point. 
Since  Kl  and  K2  must  be  positive  at  all  times,  we  account  for 
the  sign  change  in  m  and  p  in  the  different  quadrants  as  the 
calculations  proceed.  By  plotting  Kl  and  -K2  as  functions  of 
the  normalized  advance  angle,  we  can  gain  a  quick  qualitative 
feel  for  the  similarity  of  Kl  and  K2.  Figures  2.4,  2.5,  and 
2.6  are  plots  of  Kl  and  -K2  for  propellers  9,  15  and  18 
respective! y . 

As  is  readily  apparent,  even  with  a  simplified  function  for  rg 
<n)  ,  Kl  and  K2  are  y/ery  similar.  The  behavior  of  the  K<F) 
functions  in  response  to  changes  in  pitch-to-diameter  and 
expanded  area  ratios  is  also  apparent  and  are  manifest  in  the 
amplitude  variations  and  relative  smoothness  of  the  plots. 
On  the  basis  of  that  similarity,  we  form  a  new  pair  of  KO) 
functions  by  the  straightforward  transformation: 

Ka(P)  =  <K1  +  K2)/2 

Kb<P)  =  (Kl  -  K2)/2 
such   that   Kl  equals  Ka+Kb  and  K2   equals   Ka-Kb.   Hence,   Ka 
represents   an   average   of  the  two  functions  while   Kb   is   a 
measure   of  their  difference.   Figure  2.7  is  a  plot  of  the   Ka 
and  Kb  functions  for  propeller  number  15. 

A  different  approach  may  be  to  account  for  the  general  ampli- 
tude differences  in  Kl  and  K2  on  eitherside  of  the  zero  thrust 
and  torque  radials  and  themselves.   Four  gain  parameters  could 
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be  obtained  -for  each  propeller  such  that  when  Ka  and  Kb  are 
plotted,  the  maximum  amplitude  of  KbO)  will  be  considerably 
reduced.  This  was  done  -for  propeller  number  15  by  gain  scaling 
the  average  amplitudes  as  outlined  above.  Figure  2.8  shcMus  the 
result  for  Ka  and  Kb.  KbO)  now  approaches  what  may  be  con- 
sidered as  individual  propeller  nuance  characteristics. 
Since  Ka<P)  is  essentially  doubly  periodic  over  the  range  of  ^ 
we  form  yet  another  function  that  can  be  used  to  approximate 
Kl  and  K2  valid  over  the  entire  range  of  F  .  By  overlaying  the 
plots  of  KaO)  for  -135  <  P  <  45  and  45  <  F  <  225  and 
again  averaging  we  obtain  an  approximate  \«F>  function  applic- 
able to  both  the  thrust  and  torque  curves  and  valid  for  the 
two  ranges  of  P  stated  above. 

To  retain  all  of  the  information  available  such  that  we  are 
able  to  reconstruct  the  original  Kl  and  K2  functions  as 
necessary,  we  must  define  four  additional  functions.  Those 
functions  are  formed  by  first  separating  Ka  and  Kb  on  the 
intervals  of  -135  <  P  <  45  degrees  and  45  <  P  <  225  degrees. 
Labled: 

Kal  =  Ka  -135  <  P  <  45   degrees 

Ka2  =  Ka  45  <  P  <  225  degrees 

Kbl  =  Kb  -135  <  P  <  45   degrees 

Kb2  =  Kb  45  <  P  <  225  degrees. 

We  then  define  the  four  functions  that  represent  the  average 
and  differences  of  these  functions  on  the  two  intervals  of  P. 
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K      =    <Kal+Ka2)/2 

1 
K      =    <Kal-Ka2)/2 

2 
K      =:    <Kbl  +  Kb2)/2 

3 
K      =    (Kbl-Kb2)/2 

4 
Algebraic   manipulation   will    shcx«i    that: 

Kl   =   K   +K   +K   +K        -135   <    ^  <    45      deg 

12  3  4 
Kl   =   K   -K   +K   -K  45   <    P  <    225   deg 

12  3  4 
K2  =   K   +K   -K   -K        -135   <    P  <    45      deg 

12  3  4 
K2  =   K   -K   -K   +K  45    <    P  <    225   deg. 

12  3  4 
Our    approximate   KO)    -function    is   now    -found    from: 

KO)  =  <K  +K  )/2 
1   3 
valid  on  the  intervals  o-f  -135<F<45  and  45<^225  degrees. 

I-f  the  two  functions  Kl  and  K2  were  the  same  and  the  identical 

underlying  function  was  responsible  for  the  thrust  and   torque 

characteristics  on  each  side  of  the  4-quadrant  asymptote,  then 

K  ,  K  ,  and  K   would  be  zero.  In  surronary,  K<P)  is  an  average  of 

1    3        4 
Kl  and  K2  which  is  then  symmetrized  for  the  -^  and  -  torque  and 

thrust  curves.  Figure  2.9  shcxus  the  result  for   propeller 
number  9   on  the  interval  -PI  <  P  <  PI. 

To  determine  some  of  the  underlying  fundamental  characteris- 
tics of  the  K(F)  function  we  can  formulate  an  expansion  of 
K<P)  .  Two  of  the  methods  which  are  available  involve  the  Walsh 
functions  and  the  Fourier  series.  With  the  Walsh  functions, 
square  waves  are  superimposed  and  in  the  limit  can  be  used  to 
approximate  simple  sinusoids.  They  could  be  particularly  use- 
ful here  since  the  K<F)  function  does  have  ^ery  sharp  minimums 
at  the  45  and  -135  degree  angles.   The  Fourier  series  is  some- 
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what  more  familar  and  ujas  used  here.  As  a  dual  to  the  Walsh 
-functions,  in  the  limit,  the  Fouries  series  can  be  used  to 
approximate  square  waves. 

By  obtaining  the  Fourier  series  expansion  -for  KO)  we  have 
-finally,  a  potential  model  -for  this  -family  o-f  propellers. 
Table  2.3  lists  the  first  -five  Fourier  coe-f-f icien ts  for 
propellers  9,  15  and  18.  The  Fourier  expansion  is  taken  over 
the  interval  between  the  minimum  dissipation  angles.  At  those 
angles  the  K<P)  function  approaches  zero  and  provides  a 
convenient  point  to  begin  the  expansion.  The  expansion  is 
given  by  equation  <2.24>  where  c  is  defined  by  equation 
<2.25) . 

By  substituting  the  K(P)  expression  into  the  basic  modeling 
equations  (2.20)  and  (2.21),  we  can  cctfnpare  the  experimentally 
determined  thrust  and  torque  curves  to  the  model  by  plotting 
the  equations  over  the  range  of  -1  <  n  <  -♦•I  .  In  doing  so,  we 
are  assuming  that  the  differences  in  Kl  and  K2  are  sufficient- 
ly small  to  permit  their  cxnission.  Justification  of  this 
assumption  lies  in  the  accuracy  of  the  curves  generated  by  the 
model  equations  when  compared  to  actual  experimental  data. 
Figures  2.10  and  2.11  shcHfi  the  results  for  propeller  number 
18. 

KO)  =  aO  -»•  SUM<  a<k)Cos<kx)  +    b(k)Sin<kx>  )         <2.24) 
c  =  SQRT<  a3fa  -»■  b3«b  )  (2.25) 
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Figure    2. 10 
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Figure    2.11 
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Propel  1 er  Number 

9                                          15  •  18 

abc              abc  ab  c 

k=0         5.20           0         5.20         6.27           0         6.27  10.31  0  10.31 

k=l            .59      -.26         .64            .64         .59         .87  -.23  2.98  2.99 

k«2      -2.02   -1.06      2.28      -2.68   -1.15      2.92  -3.27  -1.29  3.52 

k=s3      -1.25         .27      1.28        -.57      -.19         .60  -1.25  .12  1.26 

k=4        -.59      -.87      1.06           .78      -.53         .54  -.78  .01  .78 

Table  2.3 


To  demonstrate  the  power  o-f  the  model  the  curve  resulting  -from 
utilizing  only  the  zeroth  term  o-f  the  Fourier  series  expansion 
is  plotted  along  with  the  curve  plotted  using  the  -first  three 
terms.  Doing  so  reveals  that  the  di-f-ficulty  in  -fitting  the 
proper  K  -function  lies  essentially  in  the  dissipation  regions 
o-f  the  4-quadrant  curves.  In  a  sense,  this  is  expected  since 
the  K(F>  -function  represents  a  modulating  dissipation  term  in 
the  basic  modeling  equations.  What  is  interesting  is  that  the 
plots  indicate  that  a  conwnon  dissipation  mechanism  a-f-fects 
both  the  torque  and  thrust  characteristics  and  are  lumped 
together  in  the  KO)  term. 

Additionally,  it  is  obvious  that  quite  accurate  results  -for 
normal  ahead  propeller  operation  and  -for  turbine  operation  can 
be  achieved  with  KO)  assumed  to  be  a  constant.  As  more  terms 
of  the  expansion  are  added,  the  -fit  can  be  improved  in  the 
dissipation  regions.  For  a  loss  o-f  power  simulation  -for  exam- 
ple,  the  most  accurate  results  for  thrust  and  torque  values 
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would  require  a  large  number  o-f  terms.  The  limitation  on  the 
number  o-f  terms  that  can  properly  be  added  to  the  expansion 
lies  in  the  residue  o-f  the  Kl  and  K2  functions.  Recalling  that 
KCP>  was  -formed  -from  an  averaging  process  involving  Kl  and  K2, 
the  -functions  which  represent  their  di-f-ferences  were  -formed  to 
allow  retrieval  o-f  the  original  Kl  and  K2  data.  To  the  extent 
that  these  di-fferences  degrade  the  K<P>  approximation,  their 
impact  should  properly  be  less  than  the  last  term  added  to  the 
expansion.  Again ,  the  Walsh  -functions  may  reduce  the  number  o-f 
terms  required  particularly  in  the  normal  pump  and  turbine 
areas  o-f  the  4-quadrant  diagram.  The  model  equations  can  also 
be  used  to  construct  the  -f ami  1  ar  Kt  and  Kq  versus  J  curves.  Kt 
is  given  by  equation  <2.26)  and  Kq  by  equation  <2.27) .  J  is 
simply  q/n  where  R=  n3Cn-^q3?q  and  THETA=Arctan  J. 


2 
Kt  =  l./R3{cos  (THETA)  (2.26) 


2 
Kq  =  l./R3ecos  <THETA)  <2.27) 


Here,   the   R   and  n  re-fer  to  the  appropriate  values   -for   the 
thrust  and  torque  curves. 

By  selecting  an  n  and  solving  the  modeling  equations  <2.20) 
and  <2.21)  for  q  with  a  given  KO)  ,  Kt  and  Kq  can  be  quickly 
determined  for  any  value  of  J.  Figure  2.12  is  a  plot  of  the  Kt 
and  Kq  curves  for  propeller  number  18  using  the  zeroth  term  of 
the  K<^)  expansion  and  the  first  three  terms  for  J  >  zero. 
Note  that  the  scale  factors  used  in  constructing  the   rationa- 
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lized  plots  o-f  thrust  and  torque  must  be  considered  in  calcu- 
lating Kt  and  Kq  so  that  their  ranges  o-f  values  are  those 
normally  encountered  in  the  literature.  The  scale  -factors  used 
-for  propeller  number  18  in  constructing  Firgure  2.12  can  be 
deduced  from  Table  2.2. 
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CHAPTER   3     INTRODUCTION   of   ADDITIONAL   TERMS   and   OTHER 
DETAILS  in  PROPELLER  CHARACTERISTICS 


3.  1  Pitch  as  a  primary  determinant 

As  will  be  discussed  in  Chapter  4,  the  variation  in  pitch 
expressed  as  a  pitch-to-diameter  ratio  <P/D)  has  a  significant 
e-f-fect  on  the  performance  characteristics  of  a  given  propel- 
ler. From  a  designer's  standpoint,  Saunders  C33  addresses 
pitch-to-diameter  selection  primarily  on  the  basis  of  obtain- 
ing maximum  propeller  efficiency,  particularly  when  utilizing 
series  chart  data  in  the  initial  stages  of  the  design.  From 
the  perspective  of  the  modeling  equations  however,  it  is  more 
convenient  to  see  the  relationship  between  the  relative  pitch 
to  diameter  ratio  and  specific  speed.  References  C 143  and  C 15] 
present  descriptions  on  the  use  of  the  specific  speed  para- 
meter . 

Specific  speed  variations  have  a  significant  effect  on  the  4- 
quadrant  diagram  not  easily  seen  in  the  propellers  presented 
here.  Since  propellers  are  in  general  '^ery  high  specific  speed 
machines,  the  impact  of  variations  of  that  parameter  is  some- 
what obscure.  The  principle  effect  of  specific  speed  variat- 
ions on  the  4-quadrant  diagram  was  only  touched  on  in  the 
plots  of  the  thrust  and  torque  characteristics  discussed  ear- 
lier. It  is  seen  primarily  in  the  location  of  the  zero  thrust 
and  torque  radials  of  Table  2.2.  As  the  relative  specific 
speed  increases,  the  radials  tend  to  rotate  counterclockwise 
while  the  skew  angle  addressed  in  Chapter  2  tends  to  decrease. 
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How  variations  in  pitch  a-f-fect  the  specific  characteristics  of 
a  propeller  is  best  seen  in  an  analysis  of  the  performance  of 
controllable,  reversible  pitch  propellers  <CRP) .  Reference 
1161  is  quite  thorough  in  examining  quantitatively  the  impact 
of  pitch  on  the  efficiency  and  thrust  and  torque  coefficients 
of  a  CRP  propeller.  It  is  seen  in  the  modeling  equations  as  a 
dominant  factor  in  the  determination  of  the  amplitude  of  the 
K<P)  function.  An  increase  in  the  pitch-to-diameter  ratio  has 
two  noticeable  effects.  First,  it  tends  to  increase  the  rela- 
tive amplitudes  of  the  various  harmonic  components  in  KO) 
which  results  in  a  somewhat  wider  fluctuation  of  the  4- 
quadrant  curve  in  the  principle  dissipation  regions.  Secondly, 
it  has  a  direct  affect  on  the  zeroth  term  of  the  Fourier 
expansion  which  represents  the  average  value  of  the  K<P)  fun- 
ction. The  net  effect  is  a  shrinking  of  the  4-quadrant  curves 
(corresponding  to  larger  torque  and  thrust  coefficients). 
Obviously,  the  'dc'  value  of  the  expansion  is  dominant  since 
an  examination  of  Figures  2.2  and  2.3  would  indicate  that  the 
harmonic  content  of  K<^)  in  the  dissipation  regions,  specifi- 
cally in  quadrant  I'v'  is  apparently  proportionally  equal  for 
all  three  propellers. 

3.2  Number  of  blades  and  blade  area 

The  effect  of  variations  in  numbers  of  blades  was  not  examined 
in  detail.  However,  as  pointed  out  in  Chapter  4,  the  formula- 
tion  of  the  modeling  equations  was  quite  general  and  utilizes 
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what  is  recognized  as  the  actuator  disk  theory  to  obtain  their 
initial  relationships.  A  comparison  o-f  4-quadrant  diagrams  -for 
three  bladed  and  seven  bladed  propellers  has  not  indicated  any 
potential  problems.  Mhat  is  really  the  primary  determinant  is 
the  amount  o-f  blade  area  expressed  here  as  expanded  area 
ratio.  Yet  Figures  B4(a> ,  B4(b> ,  2.2  and  2.3  would  seem  even 
to  discount  any  major  influence  o-f  EAR  on  torque  and  thrust 
characteristics    in    general. 

3.3  E-f-fects  o-f  cavitation 

Cavitation  is  probably  the  most  dif-ficult  aspect  o-f  propeller 
operation  with  which  to  deal.  Rightly  so,  a  lot  o-f  time  and 
energy  is  being  expended  to  describe  the  mechanics  of  cavita- 
tion and  its  influence  on  propeller  characteristics.  Referen- 
ces C4],  C 16] ,  £17],  CIS],  and  C 19]  explore  qualitatively  and 
quantitatively  the  role  of  cavitation  in  propeller  perform- 
ance. Reference  C 16]  relates  the  speed  of  advance  of  a  propel- 
ler to  the  cavitation  number.  By  defining  the  cavitation 
number  (sigma)  with  equation  <3.1>, 

sigma=2g<h-i-hA-hW/<v'a3^*^a  <3.  1) 

where  g  is  the  acceleration  due  to  gravity,  h  the  depth  of  the 
propeller  axis  beneath  the  water  surface,  hA  the  atmospheric 
head  pressure,  hV  the  vapor  pressure  of  water  and  Ua  the 
propeller   speed  of  advance,   it  is  shown  that  (using   typical 
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values)  Sigma  is  related  to  Va  by  the  relationship  of  equation 
<3.2) . 

si^a=  1070/VaXVa  (3.2) 

with  Va  in  Knots.  How  the  cavitation  number  a-f-fects  Kt  and  Kq 
is  seen  in  Reference  C4],  Chapter  16,  Figure  108  as  a  fall  off 
of  Kt  and  efficiency  as  the  advance  ratio  decreases.  It  does 
appear  that  the  effect  of  changes  in  sigma  for  a  higher  J  is 
smaller,  also  seen  in  Reference  C 16] . 

Yet  nowhere  is  the  role  of  cavitation  included  directly  in  the 
modeling  equations.  Since  the  intention  of  the  modeling  equa- 
tions is  to  provide  a  reasonable  method  of  determining  propel- 
ler characteristics  at  an  early  stage  of  a  ship  design  say  for 
the  purposes  of  simulation,  perhaps  an  extension  to  include 
cavitation  affects  is  unwarranted.  However,  we  can  anticipate 
where  such  affects  may  be  present. 

Equations  (2.22)  and  (2.23)  can  be  modified  to  indicate  gen- 
eral relationship  to  cavitation  number  shown  below. 

m=rg(si^a)Xq  +  Kl(P,sigma)  5^vX(n-q)  +  Rl(sigma)3(n    (3.3) 
p=rg(sigma)36n  +  K20,sigma)  3^vX(n-q)  -  R2(sigma)3«q    (3.4) 

In  some  cases,  the  presence  of  cavitation  is  handled  by  an 
effective  reduction  of  the  propeller  blade  area.  This  could  be 
incorporated  into  the  K(P)  functions.  However,  the  greater 
dissipation      would   alter    the   4-quadrant    zero    thrust    radial s    to 
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some  extent,  hence,  the  Rl  and  R2  terms  are  included  in  this 
case.  Notice  that  Rl  is  also  a  function  of  sigma  since  the 
losses  will  be  reflected  from  the  hydraulic  to  mechanial  side. 
Although  the  potential  is  there  for  inclusion  of  cavitation 
affects,  such  a  step  may  be  beyond  the  scope  of  the  assumed 
model  equations.  Again,  the  intent  of  the  equations  is  to 
provide  a  relatively  simple  means  to  simulate  propeller  per- 
formance for  first-cut  design  purposes.  Any  extension  beyond 
that  goal  Mill  increase  the  complexity  of  the  equations  rapid- 

ly. 

3.4  Hul 1 /propel  1 er  interactions 

Two  of  the  familar  parameters  which  reflect  the  influence  of 
hull/propeller  interactions  are  the  thrust  deduction  factor 
< t)  and  the  wake  fraction  <w) .  The  wake  fraction  appears  as  a 
consequence  of  the  disparity  in  ship  forward  velocity  and  the 
velocity  as  seen  by  the  propeller.  The  wake  fraction  is  de- 
fined as: 

w  =  <»v'-Va)/V 

where  Kf  is  the  ship  velocity  and  Va  the  velocity  that  the 
propeller    sees. 

The  thrust  deduction  factor  results  from  the  requirement  of 
the  propeller  to  produce  enough  thrust  to  propel  the  ship 
through   a  resistance  of  Rt  while  also  reducing  the  high  pres- 
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sure  area  over  the  ship's  stern.  De-fined  as: 

t  =  <T-Rt)/T 

where  T  is  the  net  required  thrust  to  propel  the  ship  at  a 
given  speed.  Values  for  both  -factors  are  usually  obtained  from 
model  tests.  Reference  C  16]  would  indicate  that  t  and  ui  vary 
little  over  a  wide  range  of  ship  speeds.  In  terms  of  the 
modeling  equations,  the  approach  was  one  of  considering  the 
propeller  in  an  open  water  condition  and  hence,  corrections 
for  t  and  w  were  not  included  or  necessary.  Owing  to  the 
actuator  disk  approach  essentially  taken  here,  allowances  for 
complicated  inflow  fluid  geometries  due  to  disturbances  ahead 
of  the  propeller  have  been  ignored. 

Reference  C 16]  also  addresses  the  pitch  dependency  of  the 
thrust  deduction  factor  and  gives  general  ranges  of  t  for  P/D 
ratios  of  -^  1  to  -1.  References  C4]  and  C20]  address  the  incor- 
poration of  the  thrust  deduction  and  wake  fractions  in  more 
detai 1  . 

3.5  Other  considerations 

As  formulated,  the  modeling  equatins  simply  do  not  address 
many  of  the  fine  details  of  propeller  design.  Most  of  these 
however,  including  blade  rake,  skew-back,  the  shaping  of  blade 
edges,  hub  fairing  and  blade  profiles  are  in  large  measure 
details   which   go  beyond   the  requirements   for   the   initial 

45 


stages  o-f  propeller  design.  Once  the  'macroscopic'  character- 
isics  o-f  a  generic  propeller  type  is  determined,  the  impact  o-f 
cavitation  and  speci-fic  blade  loading  requirements  can  be 
examined  but  only  in  the  detail  allcMAied  by  actual  model  tests 
or  by  a  complex  computer  program  specializing  in  propeller 
hydrodynamics. 

A  tacit  assumption  which  has  been  incorporated  -frcwn  the  begin- 
ning is  that  the  nondimensional ,  per  unit  quantity  q  is  ob- 
tainable -from  the  dimensional  values  o-f  propeller  speed  o-f 
advance  Va  and  in  turn,  the  -flow  rate  through  the  propeller. 
There  is,  hcwuever ,  a  term  called  the  axial  in-flow  -factor  a  C4] 
which  is  included  in  the  momentum  theory  -for  propeller  action 
not  directly  accounted  -for  in  the  calculations.  It  has  been 
assumed  to  be  a  quanti-f  iabl  e  constant. 

The  e-f-ficiency  o-f  a  given  propeller  and  its  relationship  to 
the  behavior  o-f  the  K<P)  function  has  not  been  explored.  The 
calculation  o-f  the  e-f-ficiency  o-f  a  given  machine  is  straight- 
forward, but  the  qualitative  relationship  of  the  efficiency  to 
pitch-to-diameter  or  expanded  area  ratio  in  the  modeling 
equations  is  yet  to  be  determined.  That  the  points  of  minimum 
dissipation  or  points  of  best  efficiency  can  be  determined 
fran  experimental  data  or  inferred  from  other  considerations 
is  accepted,  but  the  use  of  the  modeling  equations  to  estab- 
lish a  'best'  efficiency  machine  given  specific  operating 
requirements  is  a  topic  that  begs  further  work. 
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CHAPTER   4      PREDICTION  o-f  PROPELLER  PERFORMANCE   USING   THE 

MODELING  EQUATIONS 


4.  1  Comparison  with  very  similar  propellers 

Ni  th  only  a  few  o-f  the  eighteen  propellers  available  invest- 
igated in  detail  regarding  the  speci-fic  KO)  -functions,  never- 
theless it  is  possible  to  extend  the  results  to  the  remainder 
o-f  the  -family  o-f  propellers.  As  seen  by  a  comparison  o-f  pro- 
peller numbers  9,  15,  and  18,  the  KO)  -function  retains  the 
same  general  shape  over  the  range  o-f  P/D  ratios  o-f  1.0  to  1.6 
and  EAR'S  o-f  0.8  to  1.1.  Intuitively  at  least,  the  higher  EAR 
appears  to  smooth  out  the  Kl  and  K2  -functions  in  the  dissipa- 
tion regions  while  slightly  increasing  the  average  magnitude 
overall.  The  e-f-fect  o-f  EAR  is  more  easily  seen  in  the  raw  data 
plots  o-f  Appendix  B.  In  Figure  B4(b)  -for  example,  the  higher 
EAR  has  its  most  pronounced  e-f-fect  at  ^  equals  -i-PO  and  -90 
degrees.  This  is  also  shown  in  Figures  2.2,  2.3,  2.4,  and  2.5. 
A  point  by  point  comparison  will  also  indicate  that  the  varia- 
tion in  EAR  has  little  e-f-fect  on  the  thrust  characteristics 
other  than  at  those  two  particular  points  o-f  the  4-quadrant 
curve.  Seen  more  clearly  in  Figure  B4<a),  it  is  rein-forced  in 
Figures  2.2  through  2.5.  The  expanded  area  ratio  has  as  some- 
what larger  e-ffect  on  the  Kl  curve  (torque),  particularly  in 
the  normal  turbine  region  o-f  the  4-quadrant  curve.  This  is  not 
obvious  in  Figure  B4<b)  where  the  plot  o-f  the  raw  data  would 
indicate  that  EAR  has  essentially  no  e-ffect  outside  the  prin- 
ciple dissipation  regions. 
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The  e-f-fect  o-f  changes  in  pitch-to-diameter  ratio  is  much  more 
obvious.  Beginning  with  the  raw  data  o-f  Figure  B3(a),  we  can 
see  the  e-f-fect  of  P/D  ratio  throughout  the  analysis.  An  exam- 
ination o-f  Figures  2.2  and  2.3  show  this  in-fluence  quite  well. 
Note  however,  contrary  to  intuition,  the  closer  to  the  origin 
that  a  given  point  lies,  the  higher  the  relative  magnitudes  o-f 
the  corresponding  Ct  and  Cq  coe-f-f icients.  A  comparison  o-f 
Figures  2.5  and  2.6  reveal  that  the  in-fluence  of  P/D  ratio  is 
nearly  the  same  for  both  the  thrust  and  torque  characteris- 
tics. Table  2.3  points  out  that  the  various  principle  frequen- 
cies of  the  K<P)  function  change  in  amplitude  as  well  as  the 
average  'dc'  term.  The  substantial  variations  in  the  peak 
amplitudes  of  the  K<P)  function  are  also  evidenced  in  the  raw 
data  of  Appendix  B.  The  apparent  critical  points  in  Figure 
B3(a)  are  intriguing.  The  increase  in  P/D  ratio  is  reflected 
in  the  relative  specific  speeds  of  the  various  propellers  and 
is  demonstrated  in  the  gradual  counterclockwise  rotation  of 
the  zero  thrust  and  torque  radial s.  Those  apparent  critical 
points  must  lie  along  the  pivot  axis  of  that  rotation  but  the 
reason  for  the  correspondence  isn't  clear. 

By  examining  Tables  2.2  and  2.3  and  Figures  2.2  through  2.6  as 
well  as  the  figures  of  Appendix  B,  it  is  possible  to  infer 
with  \fery  reasonable  accuracy  the  appropriate  values  of  the 
rating  points  and  the  fundamental  values  of  KO)  for  all  of 
the  propellers  of  the  data  base.  With  that  knowledge,  the 
thrust  and  torque  plots  could  be  generated  and  compared  with 
Figures  2.2  and  2.3. 
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By  plotting  the  higher  harmonic  amplitudes  versus  the  zeroth 
harmonic  amplitude,  a  ccwnparison  can  be  made  between  three  o-f 
the  propellers  examined  in  detail  regarding  the  K<P)  -function, 
Figure  4.1  shows  the  amplitudes  of  the  -first  and  second 
harmonics  plotted  against  the  zeroth  harmonic  amplitude  for 
propellers  number  9,  15,  and  18.  By  extrapolating  along  the 
curves  for  a  propeller  whose  pitch-to-diameter  and  expanded 
area  ratios  would  suggest  a  particular  zeroth  harmonic 
amplitude,  we  can  approximate  the  values  of  the  Fourier  series 
expansion  terms  for  the  first  two  harmonics. 

As  an  example,  propeller  number  1  has  a  pitch-to-diameter 
ratio  of  0.6  and  an  expanded  area  ratio  of  0.5.  Assuming  the 
zeroth  harmonic  amplitude  is  equal  to  four,  the  first  two 
harmonic  amplitudes  can  be  inferred  by  extending  the  curves  to 
the  appropriate  value  of  the  zeroth  harmonic.  Using  the 
modeling  equations  with  those  amplitudes  of  the  respective 
Fourier  series  expansion  terms,  a  comparison  can  be  made 
against  the  experimentally  determined  thrust  and  torque 
curves.  Figure  4.2  shcxus  the  result  for  the  positive  thrust 
curve. 

Owing  to  the  influence  of  the  higher  harmonics  near  J  equals 
zero,  some  error  is  expected  in  that  region  of  the  thrust 
curve.  However,  as  Figure  4.2  indicates,  even  extrapolated 
values  of  the  K<^)  expansion  terms  yield  quite  reasonable 
results  for  the  normal  pump  and  turbine  portions  of  the 
curves.  The  value  of  the  zeroth  harmonic  (first  term  of  the 
Fourier   series  expansion)  has  a  substantial  influence  in   the 
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neighborhood  o-f  J  equals  zero  and  acts  with  the  neglected 
higher  order  terms  to  produce  part  of  the  error  near  the  J 
equals  zero  values  o-f  Figure  4.2.  With  more  study  o-f  the 
trends  in  the  harmonic  amplitudes  with  pitch-to-diameter  ratio 
and  expanded  area  ratio,  better  approximations  o-f  the  thrust 
and  torque  curves  o-f  similar  propellers  will  be  possible. 
Another  important  source  o-f  error  in  Figure  4.2  is  in  the  plot 
o-f  the  experimental  data.  The  scaling  which  was  done  to 
produce  the  rationalized  data  depends  on  the  values  o-f  n  and  q 
at  the  point  o-f  minimum  dissipation.  Those  scale  -factors  are 
quite  sensitive  to  errors  in  the  location  o-f  those  points, 
su-f -f icient  1  y  so  to  cause  errors  much  like  those  of  Figure  4.2. 
A  comparison  o-f  the  curves  o-f  Figure  4.2  also  points  out  a 
weakness  o-f  the  Fourier  series  representation  as  opposed  to 
other  possible  expansion  methods  in  -forming  the  K(F)  approxi- 
mation. The  absence  o-f  the  sharper  peaks  in  the  curve  gener- 
ated by  the  modeling  equations  indicates  that  a  -fairly  large 
number  o-f  terms  o-f  a  Fourier  series  expansion  is  required  to 
reproduce  the  experimental  data  in  the  -fourth  quadrant  dissi- 
pation region.  A  piecewise  continuous  approximation  o-f  K<F)  in 
the  -fourth  quadrant  would  provide  improved  results  -for  a 
simulation  requiring  better  correlation  to  experimental  data 
during  say,  transient  analysis. 
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4.2  Extrapolation  to  other  propellers 

As  developed  in  this  thesis,  the  modeling  equations  are  quite 
general  and  although  they  otAie  their  establishment  to  a  parti- 
cular set  o-f  three  b1  aded  propellers,  they  *r#  equally  Appli' 
cable  to  other  types  o-f  propellers  as  well.  Hecker  [21]  pre- 
sents 2-quadrant  data  -for  a  -few  seven  bl  aded  propellers 
including  a  hybrid.  Plotting  the  experimental  data  in  the  same 
manner  as  with  the  three  bladed  machines,  one  is  struck  with 
the  similarity  o-f  -form  o-f  the  curves  to  those  presented  earl- 
ier. Although  not  shcxun  here,  a  comparison  o-f  the  experimental 
data  to  that  generated  by  the  modeling  equations  using  only 
approximate  values  -for  the  various  parameters  was  quite  suc- 
cess-ful  .  Although  the  principle  areas  o-f  error  were  in  the 
dissipation  regions  as  expected,  the  comparison  did  not  indi- 
cate any  serious  problems  with  the  applicabilty  o-f  the  mode- 
ling equations.  A  more  detailed  comparison  o-f  the  results 
obtained  from  the  modeling  equations  to  much  di-f-ferent  propel- 
lers is  yet  to  be  done.  Considering  the  scarcity  o-f  complete 
4-quadrant  experimental  data  -for  a  large  variety  o-f  propeller 
types  and  the  limitations  inherent  in  this  -form  o-f  presenta- 
tion, such  a  ccMnparison  wasn't  possible.  However,  the  success 
o-f  the  equations  to  date  would  tend  to  en-force  the  argument 
-for  their  wide  ranging  applicability. 


53 


4.3  Extrapolation  to  missino  quadrants 

One  of  the  advantages  o-f  the  modeling  equations  lies  in  their 
ability  to  reconstruct  essentially  the  entire  4-quadrant  diag- 
ram using  only  partial  information.  Particularly  when  suffic- 
ient information  is  available  to  determine  the  point  of  mini- 
mum dissipation,  a  "universal"  K<^)  function  can  be  scaled  to 
fit  the  rationalized  data  and  the  appropriate  values  of  ^v'p  and 
ND  at  the  point  of  minimum  dissipation  substituted  to  generate 
either  the  nondimensional  Ct  and  Cq  values  or  the  dimensional 
values  of  thrust  and  torque  directly.  If  only  the  first  quad- 
rant point  of  minimum  dissipation  is  known,  the  third  quadrant 
can  be  approximated  with  little  loss  of  accuracy.  This  is 
particularly  useful  where  first  quadrant  data  exists  but  the 
model  from  which  the  data  was  taken  is  no  longer  available. 
Extending  the  concept  slightly,  a  comparison  of  available  data 
on  the  rating  points  of  various  propellers  could  form  a  data 
base  from  which  the  approximate  rating  point  of  a  specific 
machine  could  be  inferred.  Thus,  if  only  ^ery  limited  experi- 
mental results  are  available,  the  4-quadrant  curve  could  be 
constructed  with  some  confidence  by  comparing  the  modeling 
equations  results  with  whatever  information  is  available. 
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CHAPTER  5     SUMMARY 

5.  1  An  overview  o-f  the  method 

The  application  o-f  Eulerian  turbcwnachine  theory  to  the  model- 
ing and  simulation  o-f  pump/turbine  systems  has  a  long  and 
success-ful  history.  While  the  particular  method  used  here  is 
somewhat  more  recent,  this  thesis  represents  an  initial  at- 
tempt to  apply  the  methodology  to  high  specific  speed  extended 
machines.  Using  y/ery  ccwnplete  experimental  data  -for  a  -family 
o-f  three  bladed  propellers  with  varying  expanded  area  and 
pitch-to-diameter  ratios  as  a  basis,  this  procedure  was  able 
to  reveal  some  of  the  important  fundamental  behavior  under- 
lying the  dissipation  mechanisms  of  propellers  in  general. 
Perhaps  more  importantly,  the  model  equations  were  able  to 
point  out  a  common  dissipation  mechanism  for  thrust  and  tor- 
que, i.e.,  the  mechanisms  which  contribute  to  the  overall 
energy  losses  of  the  system  can  be  seen  to  ''reflect'  through 
the  transformation  of  rotational  mechanical  to  transl at ional 
fluid  energy.  The  plotting  of  the  KO)  function  for  different 
propellers  showed  their  strong  correlation  in  form,  a  correla- 
tion which  could  in  a  more  thorough  examination  be  seen  in  the 
relative  magnitudes  of  the  components  of  the  Fourier  series 
expansions.  Certainly,  the  zeroth  component  of  the  expansion 
represents  an  important  parameter  since  a  quite  reasonable  fit 
can  be  made  over  much  of  the  four  quadrants  utilizing  it 
solely  in  the  determination  of  the  model  4-quadrant  diagram. 
This   method   also  introduces  to  the  simulation   of   propeller 
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characteristics  a  procedure  which  avoids  most  o-f  the 
di-f-ficulty  with  singular  values  o-f  the  advance  coefficient.  It 
reduces  the  number  of  required  curves  to  represent  complete 
propeller  characteristics  to  only  two,  both  of  which  are 
reasonably  well  behaved  and  quite  smooth.  The  model  equations 
are  relatively  simple  and  basically  straightforward.  Thus, 
they  lend  th^nselves  to  computer  manipulation  with  little 
dif f icul ty . 

From  the  standpoint  of  the  ship  designer,  the  simulation  of 
ship  dynamics  which  include  a  realistic  propeller  in  the 
initial  stages  of  the  design  would  be  much  easier  and  certain- 
ly more  cost  effective.  A  generic  set  of  equations  with  the 
specific  values  of  a  given  propeller  rating  point  incorporated 
in  them  would  be  sufficient  to  provide  reasonable  results  in 
the  initial  design  stages  which  are  oriented  more  toward  near 
steady  state  ahead  and  astern  operation  than  a  thorough  simu- 
lation over  all  quadrants.  Once  the  geometry  of  the  propeller 
is  determined,  the  equations  can  be  updated  to  provide  in- 
creased accuracy.  That  accuracy  depends  on  the  number  of  terms 
of  the  expansion  one  is  willing  to  incorporate  into  the  simu- 
lation program.  Thus,  a  scale  model  of  the  propeller  is  not 
necessary  until  the  need  arises  for  a  self-propelled  model 
test  or  a  detailed  investigation  is  required  into  the  cavita- 
tion or  vibration  characteristics  of  the  propeller. 
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5.2  The  potential  -for  -further  applications 

Since  this  thesis  represents  an  initial  examination  oi  the 
applicability  o-f  Eulerian  turbcwnachine  theory  to  propellers, 
little  has  been  addressed  regarding  a  quantitative  ccxnparison 
o-f  various  propeller  types.  It  was  demonstrated  that  the 
results  obtained  for  the  three  bladed  propellers  which  -formed 
the  data  base  -for  this  analysis  could  be  applied  to  seven 
bladed  propellers  as  well  but  only  to  the  extent  o-f  a  quick 
qualitative  look  for  comparison  purposes. 

There  is  a  great  deal  that  can  be  done  to  enhance  the  useful- 
ness of  the  modeling  equations.  In  particular,  values  for  the 
K<P)  function  could  be  tabulated  based  on  the  examination  of 
available  propeller  data  for  various  combinations  of  governing 
parameters.  For  example,  the  first  few  terms  of  the  Fourier 
series  expansions  could  be  tabulated  for  given  pitch-to-dia- 
meter and  expanded  area  ratio  combinations  and  used  in  a  table 
look-up  scheme  such  that  a  ship  designer  wishing  to  simulate  a 
given  propeller  would  need  only  select  a  particular  rating 
point  and  combination  of  P/D  and  EAR  ratios.  The  computer 
program  could  then  select  the  necessary  coefficients  and  pro- 
ceed with  the  computations.  Or,  if  a  comparative  study  of  the 
rating  points  of  various  well  designed  propellers  used  in  ship 
propulsion  were  to  be  conducted,  a  data  base  could  be  estab- 
lished such  that  a  ship  designer  would  need  only  select  a 
generic  type  of  propeller  to  begin  a  simulation.  The  computer 
would  determine  the  appropriate  rating  point  and  corresponding 
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coe-f -f icients  to  initialize  the  program. 

The  -form  o-f  the  modeling  equations  make  them  particularly 
useful  -for  applications  to  controllable,  reversible  pitch 
propellers  <CRPP)  .  Not  only  are  they  convenient  -for  conducting 
simulations,  but  they  could  -form  the  necessary  -framework  -for  a 
control  system  synthesis.  Such  a  control  system  could  be 
use-ful  in  automated  propulsion  systems  particularly  involving 
CRP  propel  1 ers. 

There  are  other  potential  uses  including  applications 
involving  high  thrust,  low  speed  screws  like  those  used  in 
small  harbor  craft  or  tugs.  The  applications  to  exotic 
propellers  with  K^ery  high  skew  rates  or  perhaps  hybrid 
propellers  are  also  worthy  of  investigation.  Yet  another  use 
involves  an  examination  of  supercavi tat ing  propellers  and  they 
modeling  for  simulation.  With  the  promising  results  obtained 
here,  follow-on  efforts  could  be  quite  rewarding. 
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Appendix  A     Essentials  o-f  Eulerian  Turbomachines 

In  a  'black-box'  representation,  a  propeller  like  any  turbo- 
machine  can  be  modeled  as  in  Figure  Al  where  the  various 
parameters  are  torque  (t),  rate  o-f  sha-ft  rotation  (N)  ,  fluid 
velocity  <'v')  ,  and  fluid  flow  rate  <F)  .  Ne  define  an  IDEAL 
propeller  as  one  that,  in  the  process  of  power  conversion  from 
mechanical  rotation  to  axial  fluid  flow,  experiences  no  los- 
ses. Such  a  propeller  may  have  little  connection  with  reality, 
however  the  analysis  of  such  a  machine  can  provide  qualitative 
insight  into  the  mechanics  of  the  energy  transfer  process.  A 
more  comfortable  model  would  include  the  effects  of  various 
types  of  losses  into  a  PRACTICAL  propeller  model  diagrammed  as 
in  Figure  A2. 
The  losses  can  be  calculated  by  equation  Al, 

Pd  =  T^N   -  PMV  <A1) 

where  Pd  represents  the  power  dissipation.  The  amount  of 
dissipation  varies  with  the  advance  ratio  such  that  two  mini- 
mum will  occur,  one  in  the  range  of  normal  pump  operation  and 
the  other  approximately  180  degrees  away  in  the  four  quadrant 
diagram.  Using  the  positive  semi-definite  dissipation  function 
described  by  Paynter  C9],  the  equation  for  calculating  the 
power  dissipation  from  the  experimental  data  <in  nondimension- 
al  form)  is  found  to  be  simply: 

6<^)    =    2ltCqVp    -    CtJAp  <A2) 
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The  same  results  will  be  obtained  using  either  a  dimensional 
or  a  nondimensional  -formulation  o-f  the  dissipation  function. 
By  plotting  the  square  root  o-f  the  dissipation  versus  the 
angle  equal  to  the  arctangent  o-f  J,  the  quadratic  nature  o-f 
the  losses  will  result  in  a  nearly  linear  graph  such  that  the 
angle  at  which  the  minimum  dissipation  occurs  is  easily  dis- 
cernabl e. 

Figure  A3  shows  the  behavior  o-f  the  dissipation  as  a  -function 
o-f  the  advance  coe-f-f icient  -for  propeller  number  1.  (Interest- 
ingly enough,  o-f  the  two  minimum  values  o-f  "del',  the  smaller 
occurs  in  the  region  o-f  normal  turbine  operation.)  The  exten- 
sion to  a  PERFECT  Eulerian  turbomachine  comes  -from  proper 
scaling  to  eliminate  the  appropriate  amount  o-f  dissipation  so 
that  the  expression  -for  the  general  dissipation  -function  goes 
to  zero  at  the  prescribed  operating  point  (best  e-f-ficiency  or 
minimum  dissipation).  The  PRACTICAL  Eulerian  turbomachine  then 
supplements  this  dissipation  function  by  an  amount  equal  to 
the  minimum  loss. 


Per  Unit  variables  as  used  throughout  this  paper  are  defined 
as  the  nondimensional  ratio  of  a  parameter  value  to  some 
predetermined  base  value.  For  a  pump,  those  base  values  are 
taken  to  be  at  the  point  of  best  efficency  or  point  of  minimum 
dissipation.  Mhen  presenting  such  data  in  a  per  unit  form, 
obviously  the  reference  operating  point  must  be  known  to 
reconstruct  the  original  experiment. 
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Using  per  unit  variables  has  significant  advantages.  First, 
they  allcxM  free  manipulation  o-f  the  applicable  equations  with- 
out regard  to  dimensionality.  Secondly,  the  values  o-f  the  per 
unit  variables  can  be  made  to  range  in  the  neighborhood  o-f 
easily  handled  numbers  such  as  -^  1  and  -1  by  proper  selection 
o-f  the  base  values.  Thus  the  per  unit  values  give  a  qualita- 
tive as  well  as  quantitative  -feel  -for  their  behavior. 
Since  the  propeller  data  which  provided  the  basis  -for  the 
derivation  o-f  the  propeller  equations  was  already  in  a  non- 
dimensional  -form  such  that  the  physical  dimensions  o-f  Minio- 
vich's  propellers  was  no  longer  readily  apparent,  the  analysis 
was  done  in  a  non-dimensional  -form  to  retain  its  generality. 
This  may  prompt  one  to  conclude  that  the  computational  results 
are  somewhat  arbitrary,  however  one  need  only  select  'v'p,  N, 
and  D  to  derive  all  the  appropriate  values.  As  utilized  within 
the   paper,    the   per    unit    values   are   -found    by: 

p=P/Po  ,   t=T/To  ,   n=N/No  ,   q=F/Fo  ,    etc . 
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Appendix  B    Selected  Experimentally  Determined  Propeller 
Data  Converted  to  ^4-Quadrant^  Form 


The  data  presented  by  Miniovich  CIO]  and  converted  to  the 
modified  advance  coe-f-f icien t  -form  Cll]  was  converted  to  4- 
quadrant  form  by  plotting  Ct  (thrust  coefficient)  and  Cq 
(torque  coefficient)  in  polar  form  utilizing  the  convention 
shown  in  Figure  Bl.  This  is  possible  by  recognizing  that  at  a 
constant  thrust  or  torque  (which  the  Karman-Knapp  Circle  dia- 
gram utilizes),  the  defining  relations  used  by  Baker  and 
Patterson  shc^n  below: 

Ct  »  T  /  pD2(^Jp2+N2D2) 
Cq  =  Q  /  PDSCv'pS+NaDZ) 

a  =  Arc  tan  »v»p/ND 
can  be  rewritten  in  per  unit  form  as: 

R  =  1  /  SQRT(  ABS(Ct)  )    for  the  thrust  curve 


and 


R  =  1  /  SQRT(  ABS(Cq)  )    for  the  torque  curve 


and  a  =  Arctan  q/n  for  both. 

As   a   result,   the   overlap   of   the   data   in   the  ranges  of 
-1.5<J<-0.667   and   0 .667< 1/J< 1 .5   act  as  a  partial   check   by 
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plotting  on  top  o-f  each  other.  This  kind  of  representation  o-f 
the  data  yields  in  two  closed  curves  the  ccxnp-lete  characteris- 
tics of  the  propeller.  Wylie  and  Streeter  [223  address  this 
method  as  applied  to  solution  of  transient  problems  using  a 
digital  computer.  An  additional  advantage  is  the  avoidance  of 
the  problems  with  singularities  where  ever  propeller  speed  of 
advance  or  shaft  rpm  approaches  zero.  See  also  Reference  [233. 
The  curves  are  easily  recognizible  as  typical  for  high  speci- 
fic speed,  axial  flow  pumps.  See  Reference  C24]  where  Knapp 
also  describes  the  4-quadrant  representation  of  pumps  and 
turbines  and  outlines  a  procedure  by  which  the  representation 
could  be  used  to  deduce  the  transient  behavior.  Figure  B2  is  a 
sketch  of  the  flow  vs  speed  plane  indicating  the  various 
regions  of  operation.  The  zero  head  and  torque  radial s  are  the 
asymptotes  to  which  the  thrust  and  torque  curves  approach  for 
a  given  machine. 

Two  series  are  shown  here.  Figures  B3<a)  and  B3(b)  are  plots 
of  thrust  and  torque  for  propellers  1,3,  and  6.  Figures  B4<a) 
and  B4<b)  are  thrust  and  torque  curves  for  propellers  3,9  and 
15.  The  plots  reveal  the  effect  of  variation  in  pitch-to- 
diameter  ratio  and  expanded  area  ratio. 
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Figure    BSa 


67 


Figure    B3b 
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Figure    B4a 
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Figure    B4b 
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Appendix  C      The  Bond  Graph  Representation 

Paynter   C9]  presents  the  canonical  bond  graph   representation 
■for  an  Eulerian  turbomachine  as: 


GY 
RG 


Rui 


Rq 


It   is  also  represented  as  a  modulated  qyrator  series  adjoined 
to  a  shunt  resistor: 

Rs 


In   the  -first  case,   the  system  equations  are  easily  -found   to 
be: 


mssRgJfq-t-RwJ^n 


<DI) 


p=Rg3fn-Rq3Cq 


<D2> 


using    the   pump   sign    convention. 
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In    the   second   case,    the   equations   are   -found    to   be 


m  =   Rg   X    q   +    Rs   X    <n-q) 


<D3) 


p  —   Rg   3C    n    -   Rs   3f    <n-q) 


<D4) 


By   converting    <D3)    and    <D4)     to   canonical    -form,       it    is      evident 
that    Rs   can    be   expressed   as   K<P)    3C   v  where 


V  =    SORT    (    nJfn    +    qXq   ) 


<D5) 


and 


Rg  =   R6   +    Rs   3«    <n3^n    +    qXq)/v3{v 


<D6) 


A  great  deal  o-f  simplification  is  a-f-forded  with  little  loss  in 
quality  by  making  the  simple  approximation  o-f  Rg  =  n  and  Rs  = 
KJ^v  where  K  absorbs  the  individual  character  o-f  a  particular 
machine. 

As  a  more  widely  applicable  model,  general  dissipations  could 
be  adjoined  to  the  l-junctions  to  account  for  -fluid  and  mech- 
anical losses,  i.e. 


Rdn 


Rdq 
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The  equations  in  matrix  form: 


m 


Rs+Rdn 
Rg-Rs 


Rg-Rs 
Rs-Rdq 


<D7) 
<D8) 


A  more  complete  model  Mould  leave  the  dissipation  elements  as 
•functions  of  the  advance  coefficient  to  absorb  part  of  the 
individual  behavior  of  a  machine  and  smooth  out  the  require- 
ments in  the  Rs  term  to  fit  experimental  data. 

The  reader  interested  in  pursuing  the  bond  graph  representa- 
tion is  directed  to  Reference  C25]  for  an  'entry  level'  ap- 
proach to  bond  graphs  and  to  Reference  C26]  for  a  more  advan- 
ced treatment. 
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